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ELECTRON  TRAPPING  IN  SELF-MODULATED 
LASER  WAKEFIELDS  BY  RAMAN  BACKSCATTER 


Conventional  particle  acceleration  techniques  are  approaching  fundamental  limits  to  the 
accelerating  fields  due  to  material  breakdown  thresholds.  New  “structureless”  techniques  are 
now  being  investigated  to  overcome  breakdown  limitations  [1],  Some  of  these  techniques  are 
vacuum  laser  accelerators  [2],  plasma  beat-wave  accelerators  [3,4],  and  laser  wakefield 
accelerators  [3,5].  We  are  currently  investigating  the  self-modulated  laser  wakefield  accelerator 
(SM-LWFA)  concept  [6-10],  in  which  a  high  power  laser  with  a  pulse  length  much  longer  than 
the  plasma  period,  27t/o)p  -  where  C0p  is  the  plasma  frequency,  is  tightly  focused  in  a  plasma.  A 
self-modulation  instability,  caused  by  relativistic  self-focusing  (RSF)  and  the  forward  stimulated 
Raman  scattering  (SRS)  instability,  breaks  the  laser  pulse  into  beamlets  of  length  27t/(Dp  [10] 
which  resonantly  drive  a  large  amplitude  wakefield  plasma  wave.  The  plasma  wave  has  a  phase 
velocity  vp  near  the  speed  of  light,  c,  and  is  well  suited  for  high  energy  particle  acceleration.  RSF 

occurs  when  the  laser  power  exceeds  a  critical  power,  PRSF  =  17(co0/cop)  GW,  where  coo  is  the 

laser  frequency.  The  RSF  threshold  can  be  achieved  with  current  laser  technology  in  relatively 
high-density  plasmas  (n0~1019  cm'3,  PRSF~2TW).  High  density  plasmas  can  support  large 
accelerating  fields  (-100  GV/m)  before  the  onset  of  wave  breaking.  One-dimensional  (1-D)  cold 
fluid  theory  [11]  indicates  that  wave  breaking  occurs  at  an  electric  field  amplitude  of 

ewb=)/2(yp-i)Eo-  where  yp=l /^-vj/c2  and  E0  =  mco)p /e « 96-Jn 0 [cm “3 ]  V/m . 

Typically,  yp  =  G)0/cop  >  e Yp~*0  and  E0~300  GV/m  for  n0~1019  cm'3. 

Recent  experiments  have  demonstrated  electron  acceleration  in  the  SM-LWFA  regime 
[7,8],  Experiments  at  Rutherford  Appleton  Laboratory  [8]  have  achieved  background  plasma 
electron  acceleration  up  to  energies  of  100  MeV  using  a  laser  power  of  25  TW.  The  production 
of  high  energy  electrons  was  observed  to  correlate  with  a  spectral  broadening  of  the  forward  SRS 
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radiation.  This  broadening  was  attributed  to  the  destruction  of  plasma  wave  coherence  due  to 
wave  breaking.  Wave  breaking  was  used  as  the  explanation  for  the  trapping  and  acceleration  of 
background  electrons  by  the  plasma  wave  to  the  high  energies  observed. 

The  experiment  at  the  Naval  Research  Laboratory  (NRL)  was  also  conducted  in  the  SM- 
LWFA  regime.  However,  much  lower  laser  power  was  used  (2.5  TW)  and  no  evidence  of  wave 
breaking  was  apparent  although  high  energy  electrons  (up  to  30  MeV)  were  observed.  The 
experiment  used  the  laser  pulse  from  the  NRL  chirped  pulse  amplification  laser  system  focused 
in  a  supersonic  helium  gas  jet  to  generate  large  amplitude  plasma  waves.  The  1054  nm  laser 
pulse  had  a  typical  peak  power  of  2.5  TW  (400  fs  and  1  J)  and  was  focused  with  a  15  cm  focal 
length  off-axis  parabolic  mirror  to  a  vacuum  spot  radius  of  6  |im  (Ipeak  =  5xl018  W/cm2).  The  gas 
jet  used  a  supersonic  3  mm  diameter  nozzle  and  produced  a  plasma  density  of  1.4x10  cm' 
(Prsf=1-2  TW)  in  fully  ionized  helium.  The  gas  jet  density  profile  was  approximately  a  2  mm  flat 
top  with  0.5  mm  boundaries. 

An  electron  spectrometer  was  used  to  measure  the  accelerated  electron  distribution.  The 
spectrometer  consisted  of  an  electromagnet  placed  10  cm  after  the  plasma  acceleration  region 
and  a  Vi  inch  thick  plastic  scintillator  directly  coupled  to  a  photo-multiplier  tube  (PMT)  15  cm 
after  the  magnet.  The  magnet  and  scintillator/PMT  were  aligned  with  the  laser  axis.  Electrons 
with  energies  below  a  cutoff  value,  determined  by  the  magnetic  field  of  the  electromagnet,  were 
directed  away  from  the  laser  axis  and  therefore  the  scintillator/PMT.  This  “inline”  configuration 
resulted  in  measurement  of  the  integrated  number  of  electrons  above  the  cutoff  energy.  Variation 
of  the  magnetic  field  determined  the  energy  distribution  of  the  accelerated  electrons. 

The  magnet  was  a  45°  sector  magnet  with  a  2.5  mm  gap,  a  field  region  5.5  cm  long,  and  a 
maximum  field  of  2.5  kG.  Two  graphite  slabs  separated  by  a  2  mm  gap  were  placed  over  the 
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input  side  of  the  magnet  to  limit  the  acceptance  of  the  magnet  to  a  2  mm  by  2.5  mm  opening  and 
therefore  limit  the  electrons  to  a  well  defined  path  through  the  magnetic  field.  Four  inch  thick 
lead  shielding  was  placed  after  the  magnet  to  limit  the  acceptance  of  the  scintillator/PMT  to 
electrons  deflected  between  0°  and  8°.  Two  inch  thick  graphite  was  placed  over  the  lead 
shielding  to  minimize  x-ray  production  from  deflected  electrons.  The  length  of  the  magnetic 
field  region,  the  magnetic  field  strength,  and  the  acceptance  angle  of  the  scintillator/PMT  were 
used  to  calculate  the  energy  calibration  of  the  spectrometer.  This  method  results  in  an  energy 
determination  which  is  accurate  to  approximately  -5%/+25%.  The  main  source  of  error  in  the 
calculated  calibration  is  the  omission  of  the  magnet’s  fringe  fields  which  causes  up  to  a  20% 
under  estimate  of  the  electron  energy.  The  measured  electron  energies  are  therefore  the 
minimum  electron  energy  and  actual  electron  energies  are  likely  to  be  higher.  The  total  relative 
number  of  electrons  detected  at  a  variety  of  cutoff  energies  and  for  multiple  laser  shots  is  shown 
in  Fig.  1 .  The  dashed  line  represents  a  signal  to  noise  level  of  approximately  two.  Any  signal 
above  this  level  is  a  definitive  electron  detection.  Electrons  up  to  an  energy  of  28.5  MeV  were 
clearly  observed.  The  large  fluctuations  are  shot-to-shot  fluctuations  which  are  most  likely  due 
to  the  strong  non-linearity  of  the  laser-plasma  interaction  and  the  electron  trapping  mechanism. 
Measurement  of  the  continuous  electron  energy  spectrum  from  500  keV  to  5  MeV  using  film 
showed  a  smoothly  monotonically  decreasing  electron  energy  distribution  on  each  laser  shot. 

Measurement  of  the  peak  electron  energy  allows  determination  of  a  lower  bound  on  the 
electric  field  of  the  wake.  The  laser  was  observed  to  be  self-guided  through  a  plasma  over  the 
full  width  of  the  3  mm  gas  jet  [12].  The  acceleration  distance  is  therefore  the  shorter  of  the 

plasma  length  or  the  dephasing  distance  [1],  Ld  «  Yp2jtc/cop ,  which  is  the  maximum  distance  a 
particle  can  be  accelerated  before  slipping  out  of  phase  from  the  plasma  wave.  In  our 
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experiment,  Ld  ~  640  pm.  The  minimum  accelerating  field  necessary  to  generate  the  highest 
energy  electrons  observed  was  therefore  47  GeV/m. 

The  absolute  number  of  electrons  accelerated  above  1  MeV  was  also  measured.  For  this 
measurement,  the  spectrometer  cutoff  energy  was  set  to  1  MeV  and  the  scintillator/PMT  was 
replaced  with  a  silicon  surface  barrier  detector  (SBD)  with  a  detection  area  of  150  mm2.  The 
SBD  consists  of  a  500  pm  thick  disk  of  high  purity  silicon  set  between  two  metallic  plates  which 
are  biased  at  200  volts.  High  energy  particles  traveling  through  the  silicon  excited  one  electron- 
hole  pair  for  every  3.6  eV  of  energy  deposited  in  the  silicon.  Each  electron  incident  on  the  SBD 
(1-30  MeV)  passed  completely  through  the  SBD  and  deposited  approximately  400  keV  of 
energy,  i.e.,  created  l.lxlO5  electron-hole  pairs.  The  maximum  number  of  electron-hole  pairs 
created  was  measured  to  be  lxlO11.  The  number  of  electrons  passing  through  the  SBD  was 
therefore  approximately  106. 

The  number  of  electrons  measured  by  the  SBD  is  only  a  small  fraction  of  the  electrons 
accelerated  above  1  MeV  due  to  the  small  acceptance  angle  of  the  magnet.  In  order  to  determine 
the  fraction  of  electrons  entering  the  magnet,  the  beam  profile  was  measured  using  Kodak  DEF 
x-ray  film.  The  film  was  shielded  with  two  layers  of  50  pm  thick  aluminum  foil  and  a  250  pm 
thick  lead  sheet  to  block  electrons  below  1  MeV  and  stray  light.  Sixteen  shots  were  fired  to 
expose  the  film.  A  control  piece  of  film  with  only  two  layers  of  aluminum  foil  was  placed  next 
to  the  beam  profile  measurement  film,  but  out  of  the  direction  of  laser  propagation,  to  assure  that 
any  film  exposure  was  not  due  to  stray  x-rays.  No  exposure  of  the  control  piece  of  film  was 
observed.  The  beam  profile  showed  a  circular  electron  distribution  centered  on  the  laser  axis 
with  a  cone  angle  approximately  one-half  the  cone  angle  of  the  laser.  Comparing  the  electron 
beam  profile  with  the  acceptance  angle  of  the  magnet  showed  that  only  1%  of  the  electrons  were 
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detected  with  the  SBD.  The  total  number  of  electrons  accelerated  above  1  MeV  was  therefore 


approximately  10  . 

The  laser-plasma  interaction  was  investigated  by  examining  the  SRS  light  generated 
during  wakefield  production  which  consequentially  is  related  to  the  wakefield  amplitude.  For 
large  amplitude  plasma  waves,  non-linear  steepening  results  in  harmonics  (coo  ±  ntOp)  of  the 
standard  SRS  spectrum  (coo  ±  CDp)  [13].  SRS  light  was  examined  by  placing  a  lens  with  an 
acceptance  angle  of  ±10°  at  0°,  20°,  30°,  40°,  60°,  and  90°  to  the  laser  axis  to  image  anti-Stokes 
light  into  a  0.25  m  spectrometer.  The  spectrometer  measured  the  scattered  spectrum  from 
approximately  700  nm  to  1000  nm  which  allowed  the  first  four  harmonics  of  the  anti-Stokes  line 
to  be  measured  on  each  laser  shot.  The  anti-Stokes  spectra  from  0-60°  were  qualitatively  similar 
with  multiple  orders  and  no  broadening  of  the  anti-Stokes  peaks  apparent.  No  anti-Stokes  signal 
was  observed  at  90°.  The  angle  of  40°  was  found  empirically  to  provide  the  highest  signal-to- 
noise  ratio  measurement  of  the  SRS  light  and  therefore  the  wakefield  amplitude.  Fig.  2  shows  a 
typical  spectrum  at  40°  where  up  to  the  4th  harmonic  is  visible.  The  5th  harmonic  was  also 
measured  on  many  shots  by  lowering  the  central  wavelength  of  the  spectrometer.  The  large 
angular  spread  of  the  anti-Stokes  light  is  most  likely  due  to  the  highly  3-dimensional  nature  of 
the  wakefield  due  to  the  tight  focusing  of  the  laser  pulse  (focal  spot  size  ~  plasma  wavelength) 
[14,15].  The  multiple  orders  of  SRS  observed  and  their  comparable  intensities  indicates  that  a 
highly  non-linear,  large  amplitude  plasma  wave  was  present.  The  non-linearity  of  the  plasma 
wave  was  correlated  with  the  high  energy  electron  production  (see  Fig.  2  inset),  where  the 
intensity  of  the  2nd  harmonic  exhibited  a  strong  correlation  with  the  number  of  accelerated 
electrons.  In  no  instance  during  the  experiment  was  any  broadening  of  the  scattered  peaks 
observed.  For  example,  the  width  of  the  1st  harmonic  shown  in  Fig.  2,  when  the  plasma  wave 
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was  highly  non-linear,  remained  the  same  at  lower  laser  powers  when  the  plasma  wave  remained 
linear.  This  indicates  that  the  plasma  waves  were  highly  coherent  with  no  evidence  of  wave 
breaking.  Large  shot-to-shot  fluctuations  were  observed  in  the  SRS  spectrum  concurrent  with 
energetic  electron  production.  The  SRS  spectrum  and  electron  production  fluctuations  did  not 
correlate  with  laser  power,  gas  jet  density,  or  any  other  directly  controllable  experimental 
parameter,  except  that  for  very  low  laser  power  (  <  1  TW  )  or  low  plasma  density  (  <  10  cm'  ), 
high  energy  electrons  were  never  produced. 

The  generation  of  high  energy  electrons  in  this  experiment  without  evidence  for  wave 
breaking  strongly  suggests  that  electrons  are  pre-accelerated  prior  to  being  trapped  by  the 
wakefield.  Previous  experiments  on  our  system  have  shown  that  approximately  10%  of  the 
incident  laser  pulse  is  reflected  by  backward  Raman  scattering  (BRS)  [16].  One  candidate  for 
pre-accelerating  electrons  is  the  low  phase  velocity  waves  generated  in  BRS  [17].  BRS  produces 
backward  traveling  light  of  frequency  ®o  -  A®,  where  ®o  »  A®  >  ®p.  The  BRS  light  wave  can 
beat  with  the  main  laser  pulse  to  generate  a  low  phase  velocity,  vpb=c(3pb  »  (Acq/2®0)c, 

forward  traveling  beat  wave.  At  high  laser  intensities,  BRS  occurs  in  the  strong-pump  limit,  and 
the  effects  of  the  space-charge  wave  can  be  neglected  [18].  The  low  phase  velocity  beat  wave 
can  pick  up  some  of  the  low  energy  electrons  in  the  background  plasma  and  accelerate  them  to 
sufficient  energies  so  as  to  be  trapped  by  the  high  phase  velocity  wakefield  (vp  =  c).  These 

trapped  electrons  can  then  be  accelerated  to  much  higher  energies  by  the  wakefield. 

A  1-D  numerical  simulation  of  this  two  stage  acceleration  mechanism  shows  the 
qualitative  behavior  observed  in  the  experiment.  The  simulation  prescribes  three  analytical  wave 
potentials  (normalized  to  mc2/e)  corresponding  to  the  primary  laser  pulse  ao,  the  wakefield 
plasma  wave  <)),  and  the  BRS  electromagnetic  wave  a,,  and  then  pushes  test  particles  in  these 
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potentials.  The  simulation  uses  the  coordinates  C,  =  z  -  ct  and  x  =  t,  the  potentials  are  initialized 
in  the  region  C,  <  0  (C,  =  0  corresponds  to  the  initial  front  of  the  laser  pulse),  and  the  particles  are 
initialized  at  rest  in  the  region  C,  >  0. 

The  laser  pulse  parameters  used  in  the  simulation  were  chosen  to  closely  model  the 
experimental  laser  pulse  parameters  —  peak  normalized  amplitude  aom  =  1.4,  wavelength 
A,  =  1  pm,  and  normalized  frequency  C£>o/C0p  =  8.5.  The  BRS  is  assumed  to  saturate  at  a 
normalized  amplitude  of  aim  =  0.033,  based  on  analytical  estimates  for  saturation  in  the  strong- 
pump  regime  [18].  All  potentials  have  a  characteristic  rise  length  of  Tr  =  25/(Op .  Fig.  3  plots  the 

normalized  axial  momentum,  uz  =  ypz  =  pz/mc,  of  simulation  electrons  at  a  wakefield  potential  of 
<t)o  =  0.6  as  a  function  of  the  coordinate  t,  after  250  pm  of  propagation.  This  demonstrates  the 
initial  trapping  of  particles  in  the  combined  wake  and  BRS  fields. 

Production  of  energetic  electrons  is  seen  in  the  simulations  only  when  the  wakefield  <t>o 
and  BRS  ai  amplitudes  are  sufficiently  large.  Fig.  4  plots  the  peak  electron  energy,  Wmax,  and 
percentage,  ftr,  of  plasma  electrons  which  are  trapped  and  accelerated  to  energies  exceeding  10 
MeV  as  a  function  of  <J>o  after  propagating  2.5  mm  (other  parameters  are  as  in  Fig.  3).  The 
threshold  for  trapping  is  at  (j)0  =  0.55,  beyond  which  ftr  rapidly  increases.  Once  trapping  occurs, 
Wmax  corresponds  to  the  detuning  limit,  i.e.,  Wmax  =  4y2<|)0mc2.  This  is  in  qualitative  agreement 

with  the  strong  correlation  of  wakefield  amplitude  to  high  energy  electron  production  observed  in 
the  experiments  (see  Fig.  2  inset). 

An  estimate  for  the  trapping  threshold  can  be  obtained  by  requiring  that  the  wakefield 
separatrix  overlap  the  beat  wave  separatrix  [1,19].  In  momentum  phase  space  (uz  versus  z-vpt), 

the  minimum  of  the  wakefield  separatrix  is  given  by  uw  min  =  y2  /4§0  -<|)0  ,  assuming 
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2<t>oYp/lL  » 1  where  y±  =  yi  +  aQ.  Likewise,  the  maximum  of  the  beat  wave  separatrix  is 
given  by  ub  max  «  ppbYj_  +2>/a0a] .  Passing  of  electrons  from  the  slow  beat  wave  to  the  fast 

wakefield  can  occur  when  Ub.max  ^  uw,mjn,  i.e.,  <J)0  ^(yjl  _ub,max)/2-  This  condition  gives 

4>o  >  0.6  for  the  parameters  of  Fig.  3  (f3pb  =  0.059,  ao  =  1.4,  and  ai  =  0.033).  A  more  accurate 
calculation  gives  a  trapping  threshold  of  <J)0  >  0.54,  in  excellent  agreement  with  the  simulation 
results.  Hence,  self-trapping  of  plasma  electrons  can  occur  when  (po  <  1,  i  e.,  well  below  the 
theoretical  wave  breaking  amplitude  of  <j)0  =  EWB/E0  *  4. 

In  conclusion,  we  have  observed  very  high  energy  electrons  (up  to  30  MeV)  accelerated 
in  a  SM-LWFA  at  relatively  low  laser  power  (2.5  TW).  Optical  diagnostics  show  the  existence 
of  highly  nonlinear,  large  amplitude  plasma  waves.  Using  the  dephasing  distance  as  the 
acceleration  length,  the  acceleration  gradient  is  estimated  to  be  greater  than  47  GeV/m.  The  high 
energy  electrons  are  observed  with  no  sign  of  wave  breaking.  Numerical  simulations  show  that 
low  phase  velocity  beat  waves  generated  by  BRS  can  pre-accelerate  background  plasma  electrons 
to  sufficient  energy  to  be  trapped  by  the  high  phase  velocity  wakefield.  Both  the  experiment  and 
the  simulation  show  that  the  trapping  and  acceleration  of  electrons  in  the  SM-LWFA  is  a  highly 
non-linear  process  strongly  dependent  on  the  wakefield  amplitude. 
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Fig.  1  :  Electron  energy  distribution  measured  with  the  scintillator  and  PMT.  The  dashed  line 
represents  a  signal  to  noise  level  of  approximately  two.  Any  signal  above  this  level  is  a  clearly 
discernible  electron  peak. 
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Fig.  2  :  The  multiple  order  anti-Stokes  spectrum.  The  well  defined  peaks  show  a  highly 
nonlinear  and  coherent  plasma  wave  with  no  evidence  of  wave  breaking.  The  inset  shows  the 
relationship  between  the  2nd  harmonic  intensity  and  the  number  of  electrons  accelerated  above  1 
MeV. 
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Fig.  3  :  Momentum  phase  space  plot  of  the  test  electrons  in  the  simulation  after  propagating 
ct=250  pm  with  ao=1.4,  ai=0.033,  <|>o=0.6,  and  (Oo/(Op=8.5.  The  laser  pulse  resides  in  the  region 
z  -  ct  <  0  and  propagates  to  the  right. 
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Fig.  4  :  Simulation  results  showing  the  maximum  energy  (solid  line)  and  fraction  above  10  MeV 
(dashed  line)  of  accelerated  electrons  as  a  function  of  the  wakefield  potential  0o  at  cx=2.5  mm  for 
the  parameters  of  Fig.  3.  The  y-axis  is  in  units  of  MeV  for  the  energy  and  the  percentage  for  the 
fraction. 
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